Theory of Electron Spin Relaxation in ZnO 
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Doped ZnO is a promising material for spintronics applications. For such applications, it is important to 
understand the spin dynamics and particularly the spin relaxation times of this II-VI semiconductor. The 
spin relaxation time x s has been measured by optical orientation experiments, and it shows a surprising non- 
monotonic behavior with temperature. We explain this behavior by invoking spin exchange between localized 
and extended states. Interestingly, the effects of spin-orbit coupling are by no means negligible, in spite of 
the relatively small valence band splitting. This is due to the wurtzite crystal structure of ZnO. Detailed 
analysis allows us to characterize the impurity binding energies and densities, showing that optical orientation 
experiments can be used as a characterization tool for semiconductor samples. 



I. INTRODUCTION 

Zinc Oxide has been the subject of considerable experimen- 
tal and theoretical investigation for many years. 1 Its bandgap 
is in the near ultraviolet, making it useful as a transparent con- 
ductor and as sunscreen. Its piezoelectricity opens up trans- 
duction applications. The activity has intensified more re- 
cently because of the possibility that ZnO might be useful for 
spintronics or spin-based quantum computation. It has been 
predicted to be a room-temperature ferromagnet when doped 
with Mn. 2 Furthermore, its spin-orbit coupling is generally 
thought to be very weak compared with GaAs. The usual 
measure of the strength of spin-orbit coupling in semiconduc- 
tors is the energy splitting at the top of the valence band. It 
is said that the spin-orbit coupling is negligible in ZnO be- 
cause the valence-band splitting is —3.5 meV, 3 as opposed to 
340 meV for GaAs. Smaller spin-orbit coupling should lead 
to long spin relaxation times. Long relaxation times are re- 
quired if spin information is to be transported over appreciable 
distances. 

The spin relaxation time x s has been measured by Ghosh 
et al4> to be about 20 ns from to 20 K in optical orientation 
experiments. X s is sometimes called T 2 * even in the absence 
of an external field. Since the data from Ref. (4) used in this 
paper were taken at zero field, the relaxation time is taken to 
be x s to avoid confusion with experiments conducted at finite 
field. The data show two surprising features. First, the relax- 
ation times are actually somewhat shorter than the longest re- 
laxation times in GaAs, which are about 100 ns. 5 One might 
expect the opposite given the relative strength of spin-orbit 
coupling in the two materials. Second, x s shows a non- 
monotonic temperature dependence, first increasing slightly 
and then rapidly decreasing - but increasing temperature usu- 
ally promotes spin relaxation. 

In this paper, we show that a theory previously developed 
for T s in GaAs 6 can account for these observations. The the- 
ory must be modified to take account of the different impurity 
levels and binding energies of ZnO. This is important, be- 
cause, in spite of intensive investigation, the nature of the im- 
purities that govern the electrical properties of ZnO remains 
controversial, and our analysis sheds some light on this issue. 
Even more interestingly, it turns out that the wurtzite crystal 



structure has very important consequences for the D'yakonov- 
Perel (DP) 7 scattering that dominates the relaxation at higher 
temperatures. Thus the crystal structure must be taken into 
account fully. The final message will be that the "weak" spin- 
orbit coupling of ZnO is not negligible for spin relaxation, and 
it does not lead to long relaxation times. 

In the next section we give the background information for 
ZnO. Sec. Ill is devoted to a derivation of the equations of 
motion for the spins. In Sec. IV the computational method is 
described. The determination of the parameters in the equa- 
tions of motion is a separate task. The most important of 
the parameters is that which controls that DP spin relaxation. 
Since the calculation of this parameters is not straightforward, 
we devote Sec. V to that. Sec. VI gives the results. Sec. VII 
is the conclusions and puts the results into context. 



II. BACKGROUND 

In ZnO produced by the hydrothermal method, it is gen- 
erally thought that there are two sets of impurity states, 
one shallow and quasi-hydrogenic, one deep and very well 
localized^ Their precise physical nature is not known. In 
the case of the deep impurity, it is believed that a lattice defect 
accompanies the chemical impurity. The binding energies are 
in the range of a few 10s of meV for the shallow impurity and 
a few 100s of meV for the deep impurity. We shall demon- 
strate below that the optical orientation data can put bounds 
on these numbers. 

ZnO crystallizes in the wurtzite structure rather than the 
zincblende structure familiar from the III-V compounds. This 
has very important implications for the conduction band 
states. The spin-orbit interaction lifts the spin degeneracy 
in the conduction band. In zincblende structures crystal sym- 
metry implies that the splitting is cubic in the magnitude of 
the wavevector k, but in the wurtzite structure the splitting is 
linear. 10 However, the spin relaxation time of the low-lying 
conduction band states depends mainly on the spin splitting 
near the conduction band minimum, and this is larger in ZnO 
than in GaAs for small enough k. 

In optical orientation experiments, electrons are excited 
from the valence band to the conduction band by circularly 
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polarized light tuned close to the bandgap energy (pump 
step). The population of conduction electrons so created is 
spin-polarized.— Energy relaxation then occurs on a short 
time scale (< 1 ns), but most of this relaxation is from spin- 
conserving processes, so there is a longer time scale (or time 
scales) on which the spin of the system relaxes. When an 
external magnetic field is present, the time scale to relax the 
transverse component of the net magnetization is called x s . It 
is measured using Faraday rotation or the Kerr effect (probe 
step). 

The important physical point is that the fast energy relax- 
ation leads to a thermal charge distribution for the electrons 
by the time 1 ns has elapsed, but the spin distribution relaxes 
on longer time scales. The thermal charge distribution means 
that the localized donor impurity states are mostly full at the 
relatively low temperatures of the experiment. The spins of 
the localized electrons must be included along with the con- 
duction electron spins. The spins of localized and extended 
states can be interchanged by the exchange coupling, a pro- 
cess we call cross- relaxation. This is often a rather fast pro- 
cess and is particularly important when the relaxation times 
of the localized and extended states are very different in mag- 
nitude. In GaAs this process is important in all the regimes 
of temperature, applied field, and impurity density that have 
been studied, and it is important in ZnO as well. 

In the following section we derive a set of modified Bloch 
equations to describe the aforementioned spin dynamics. 



III. MODIFIED BLOCH EQUATIONS 



~ og of the impurity. When |r — R/| > as, the conduction 
electron proceeds unhindered by the effective field. This ef- 
fective field is a result of the exchange potential. An itinerant 
electron will spend an average time of ag/v within the range 
of the effective field where v is the velocity of the electron. 
Thus the time between encounters 1 ^ is 1 /«/a|v. 

In a semiclassical picture the spin of the itinerant electron 
undergoes precession of magnitude A(j) = Jas/2v through a 
random angle during each encounter with an impurity. The 
spin of the impurity electron also precesses but with angle 
— A(|). Since the sum of spins, s c + S/, commutes with +H[, 
the total spin in the system must be conserved. However the 
spin in each subsystem may shift between one another; this is 
cross-relaxation. 

It turns out for the parameters of the system under consid- 
eration that A(j) ~ 1 , and we then find that 



1 



n,a|v 



(2) 



which implies that the spin is essentially randomized after one 
impurity encounter. 

If we consider an ensemble of conduction electrons with 
a net magnetization m c , this magnetization is exchanged at 
a rate of As previously mentioned, any magnetiza- 

tion lost from the conduction electrons must be gained by 
the localized electrons and vice-versa. For clarity we write 
1 j% c J = ni/Y r and 1 /if — n c /Y'' where f ' r = 1 /a 2 B v. 

We now examine the second term of the Hamiltonian 



g/j B (b x (t)a x + b y (t)a y + b z (t) a z ) 



(3) 



We consider a conduction electron in the semiclassical ap- 
proximation. It moves as a wavepacket with a well-defined 
momentum and scatters from impurities and phonons at time 
intervals of average length z p , where z p is the momentum re- 
laxation time. Its spin operator is s c . The spin-dependent part 
of its Hamiltonian in the absence of an external magnetic field 
is: 

//'■ = Hf + HI = - JL £ / ( r _ R,.) S/ . Sc _ g B h {t) . Sc . (l) 

The first term, ffj, is the exchange interaction with impurity 
spins s, located at an positions R,. It is the same interaction 
that is responsible for the Kondo effect, but the temperatures 
here are all much greater than the Kondo temperature. The 
range of the function J (r — R,) is roughly as, where as is the 
effective Bohr radius. The second term, HS, represents other 
spin relaxation mechanisms that we model as a small random 
classical field b(f ) with a correlation time much shorter than 
Z s . An analogous Hamiltonian H 1 can be written for a local- 
ized electron. 

First, we concentrate on the spin dynamics resulting from 
the spin-spin term and ignore the second term. In the dilute 

1 /3 

limit iflBfii' <C 1), a conduction electron encounters impu- 
rities with randomly aligned spins if no short-range order is 
present in the impurity system. An effective field from the 
impurity spin affects the conduction electron when it is within 



This Hamiltonian relaxes the conduction electron spin. To 
extract a relaxation rate from this Hamiltonian, we use the 
equation of motion 



dp(t) 
dt 



■Mt),Hm] 



(4) 



where p(f) is the 2x2 spin density matrix for an electron of 
a given momentum. We assume that the total density matrix 
for the conduction electron factorizes; we neglect off-diagonal 
terms that come from correlations. By iteration, we can write 
this equation as 

dp jt) 
dt 

^<[p (0) M (01) - £ /„' <[[P (0 M (0] ' H 2 (')]>*' 



(5) 



where the angular brackets indicate averaging over all ori- 
entations of b(t). To simplify notation, from now on angu- 
lar brackets will be suppressed on the density matrix. Since 
(bi(t)) = 0, the first term is zero. We assume that differ- 
ent directions of bi are uncorrected and (since the external 
field is zero) different direction are equivalent. Then we have 
{bi{t)bj{t')) = (b(t)b(t'))8ij. Therefore, Eq. © reduces to 



dp (t) 
dt 



2 2 
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Yj_p{t'),Vi}*i{b(t)b{t'))dt'. (6) 
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The correlation function is assumed to be stationary in time so 
(b(t)b(t')) = g(t' —t)= g{l)£- If the correlation time of the 
^-fluctuations, i e , is short, p will not change on that timescale 
and g(z) will be nearly a 8- function (Markov approximation). 
Eq. © can then be written as 



dp (t) _ 2gY B 1 



dt 



H 2 4 



z £[p(0,°/]<*y o (b(t)b(t>))dt>. (7) 



The integral is approximated by (b 
time scale T c by 

X c \ h J 



)i P . Define the relaxation 



(b 2 )Xe 



giving 



dp(t) 1 ^ 



(8) 



(9) 



The density matrix can be expanded in Pauli spin matrices 

(10) 



p{t) = ^I+W' n '( t ) a i- 



where / is the 2x2 identity matrix and m; = 7>(G,p) is the 
expected value of the magnetization. Inserting Eq. ( TTOb in 
Eq. (O and matching coefficients of Pauli matrices gives 
a set of equations for the dynamics of m. For instance 
for conduction electron magnetization m c in the x-direction, 
dm c /dt = Tr(a x dp/dt) = —m c /x c . As with Hf, similar ex- 
pressions for the localized magnetization m/ can be found: 
dmi/dt = Tr(o x dp/dt) = — m//T/. 

By combining the effects of Hi = H[+H[ and H2 = #| + 
Hi, the modified Bloch equations for the magnetizations can 
be expressed as 



ni_ 

r 



dm c f 1 

dt \%c 

dm) n\ 

= — m c — I — 

dt f r V T; 



n c 

i c H mi 

r) m <- 



(id 



for two spin systems - itinerant and localized spins. T c and 
T; in Eq. (fTTT i are now described in terms of well known re- 
laxation mechanisms which will be discussed in the next sec- 
tion. This model was successfully applied to GaAs. 6 For ZnO, 
these Bloch equations are easily extended to account for the 
multiple-type impurities present. 



IV. METHOD 

We now seek to write equations like those of Eq. ( fTTT ) with 
regard given to the two types of impurities in ZnO - shal- 
low and deep. As mentioned above, we find that the cross- 
relaxation is important to understand the data. These rates 
come from the Kondo-like Jsi ■ s c interaction between an im- 
purity spin S/ and a conduction band spin s £ . An expression 
for / in terms of tight-binding parameters can be derived us- 
ing the Schrieffer-Wolf transformation. 14 One expects that 



the cross-relaxation between conduction and shallow donor 
electrons to be much more rapid than the cross-relaxation be- 
tween conduction and deep donor electrons because of the 
greater binding energy of the deep impurity and its larger on- 
site Coulomb energy. This is confirmed by the fit to the data. 
In fact we find that terms involving cross-relaxation between 
the deep donors and either the conduction band electrons or 
the shallow donor electrons can be neglected. With these sim- 
plifications, for ZnO Eq. (fTTT i extends to 



dm c ( 1 «/,. \ n c 
= — m r H m/. 

dm u n h / 1 n c \ 



dt r c r iS 

dmu 1 



(12) 



dt 



-mid- 



In this equation, ra c , m; s , and m^ stand for the magnetizations 
of the conduction electrons, the electrons on shallow impuri- 
ties, and the electrons on deep impurities, respectively. The 
n's denote the corresponding volume densities. Each of the 
populations has a relaxation time T C ,T/ S , and X«. From Eq. 
(fT2l . we can then find the magnetization as a function of time. 

Standard methods can be used to solve these differential 
equations. The solutions yield a time dependence of the total 
magnetization, m(t) = m c (t) + mi s (t) + »J/ c /(f )> to be a sum of 
three exponentials, exp(— T + t), exp(— T_f), and exp(— Tjt) 
where 



(13) 




with S given by 



\ 




4n c n /j 
7S 2 



(14) 



No net moment can exist on the deep donor sites since no 
moment is excited into the deep states on account of them 
being significantly below the conduction band, and no net 
moment cross relaxes into these states. Therefore Y c ( can be 
ruled out as being the observed relaxation rate. In the regime 
that (n/j + n c )/^ r s ^> 1/t c , 1/t/ 4 , the rate F + simplifies to 
(n c + n/i)/7^ T s and is very rapid and the rate F is slower, 



1 



1 



n c + n h % c n c + n h x u 



(15) 



We fit the data with this equation and associate it with x s . 
We see that the relaxation rate depends on two factors: the 
thermodynamic occupations of the shallow donors (the deep 
donors are always nearly full in the temperature range studied 
here) and form of the relaxation rates for the conduction and 
localized shallow states. 

The densities can be computed using standard formulas 
from equilibrium statistical mechanics, since we deal only 
with time scales long compared to the fast energy relaxation 
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scale. As a function of temperature T, the ratio n c /ni s natu- 
rally increases rapidly as T — > \£i s \/ks, where £/ f is the bind- 
ing energy of the shallow impurity. \zid\ is so large that these 
states are always occupied at the experimental temperatures, 
which range from 5K to 80K. 

T c is fairly complicated to calculate because there are 
several mechanisms that can relax the conduction electron 
spins. The simplest such mechanism is the Elliot- Yafet (EY) 
process 15 that arises from spin mixing in the wavefunctions. 
When a conduction electron is scattered by a spin-independent 
potential from state k to state k', the initial and final states 
are not eigenstates of the spin projection operator S z so the 
process relaxes the spin. The rate of relaxation due to the 
EY process is well known to be of the form: 1/xey — 
V-eyT 2 /x p (T) where <Xey is a material dependent parame- 
ter and T p is the momentum relaxation time. 16 We estimate 
a E y(th) = 4.6 x 1(T 15 K~ 2 . The Bir-Aronov-Pikus (BAP) 
mechanism arises from the scattering of electron and holes. 
This relaxation mechanism is commonly considered to be 
negligible in n-type materials like those under consideration 
here since the number of holes is small. 18 The D'yakonov- 
Perel' (DP) mechanism 7 arises from the ordinary scattering 
of conduction-band states. Since this has not previously been 
calculated in a wurtzite structure, we devote the next section 
to it. This calculation yields an expression for i c as a function 
of temperature. 

T/ s and Zid are due to non-spin-conserving anisotropic 
exchange (Dzyaloshinski-Moriya) interactions ! 19 ' 2Q The 
anisotropic exchange term is important. It arises from 
spin-orbit coupling and produces a term proportional to 
d • Si x S2 where d is related to the interspin separation and the 
exchange integral between the wavefunction on sites 1 and 
2. However, it is not possible to calculate it in detail when 
the nature of the impurities is not well known. We estimate 
the rate as \/i DM = a DM (ni mp , s + n impA ) where n imp . s and 
nim P ,d are the total impurity concentrations of the shallow and 
deep impurity respectively and (Xdm has a weak temperature 
dependence that we neglect. The main contribution comes 
from the the overlap of the shallow impurity wavefunctions, 
which we take to be hydrogenic, with the deep impurity 
wavefunctions, which we take to be well-localized on an 
atomic scale. The details of how to estimate the resulting 
relaxation may be found in Refs. (6.211221). The numerical 
value we find from theory is aoM{th) = 1.12 x 10~ 20 cm 3 
ns~'. When nuclei possess nonzero magnetic moments, the 
hyperfine interaction between electron and nuclear spin is a 
source of spin relaxation for localized electrons. 23 However, 
zero nuclear spin isotopes of Zn and O are 96% and 99.5% 
naturally abundant respectively. Therefore we rule out the 
hyperfine interaction from being an observed relaxation 
mechanism in Ref. (0). 



V. DP MECHANISM IN WURTZITE CRYSTAL 
STRUCTURES 

The conduction band states undergo ordinary impurity and 
phonon scattering. Each scattering event give a change in the 



wavevector k, which in turn changes the effective magnetic 
field on the spin that comes from spin-orbit coupling. This 
fluctuating field relaxes the spin. The effective field strength 
is proportional to the conduction band spin splitting. Bulk 
zincblende crystals have conduction band splittings cubic-in-fc 
due to bulk inversion asymmetry ( Dresselhaus effect )i2i In 
addition to cubic terms, bulk wurtzite conduction bands also 
possess spin splittings proportional to linear terms in k due 
to the hexagonal c axis which gives bulk wurtzite a reflec- 
tion asymmetry similar to the Rashba effecti 10 - 25 i 26 - 27 We can 
write the spin-orbit Hamiltonian to include both the Rashba 
and Dresselhaus terms: 

ff io (k) = [<XKi + YK 3 ]-a (16) 

where Ki = (k y , — k x , 0) is linear-in-£, K3 = (k 2 — 

bk 2 ){k y , — k x ,0) is cubic-in-fc, a = (o x ,a y ,o z ) are the Pauli 
spin matrices, and a, y are spin splitting coefficients. 1Q i 26 i 28 
The parameter b is roughly equal to four for all wurtzite 
materials^ Note that there is no spin splitting along the 
hexagonal axis (z). 

The linear-in-fc term dominates and we can determine the 
spin relaxation rate by following the treatment given by Pikus 
and Titkov in Ref. (I29l) which yields the following relaxation 
rates: 

1 4a 2 

X DP,ii 



1 4a 2 

where the overbar denotes angular averaging, and i,j denote 
the Cartesian components of Ki. The momentum relaxation 
rate is defined as 

1 r l 

— =/ o(0)(l -cos9)<icose (19) 
t P J-i 

where 0(8) is the scattering cross section and is the an- 
gle between initial and final ki 2 ^ In bulk wurtzite Kijj 2 = 
ici^ 2 = Kl7 2 = k 2 /3, K1T 2 = 0, and in the unstrained crystal, 
KijK\j = for ; 7^ /'. From Eq. ( fTTb , we can write 

1 _ 2 _ 4 a 2 2 _ 8 OT *« 2 - F 

JTT - JIT - 3tf z p k -3T Tp k (20) 

L DP,|| l DP,z 

where m* is the electron effective mass and is the energy 
h 2 k 2 /2m*. This result can be Boltzmann averaged (denoted 
by angle brackets) to obtain 

^=feH^ £i) = a ^ (r,r <21) 

where a£] = 4-m*a 2 k B /h 4 and z p (T) = (f p E k )/(E k ) = 
2{z p Ek) /liksT . The temperature dependent momentum re- 
laxation time, 1 P {T), can be determined from electron mobil- 
ity (ji e ) measurements from /j e = ex p (T)/m* where e is the 
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charge of an electron, a has been calculated 10 to be 1.1 x 

1(T 4 eV-nm which gives a theoretical value of ttppfyh) — 34.6 
KT 1 ns~ 2 . 

Similarly, the cubic-in-£ term can be calculated to be 



1 



1 



i (3) 

l DP 



(T) 



DP, 



l DP, 



= a^x p (r)r 3 (22) 



where a^p = E0Qj 2 m* 3 k p /3h & where the dimensionless 
quantity Q depends on the type of scattering and is of order 
unity, yhas been calculated 10 to be 3.3 x 10~ 4 eV-nm 3 which 

yields a$ P {th) = 2.0 x 10~ 4 K~ 3 ns~ 2 . 

The sample from which the momentum relaxation times 
l p (T) were extracted^ was hydrothermally grown by the 
same company as Ghosh et al. sample in Ref. dj). 



VI. RESULTS AND DISCUSSION 

In Fig. 1 we show that temperature dependence of x s as mea- 
sured in a bulk ZnO sample and our fit (using Eq. (JT3J) to 
the data. It is seen immediately that the temperature de- 
pendence is not monotonic and that this is well-reproduced 
by the theory. The reason is simple. At low temperatures 
T <C |£/. s |/&b nearly all the electrons are in localized states. 
These states relax by the temperature-independent DM mech- 
anism: 1 /lis = 1/tdm- This mechanism alone determines the 
T = values. When T approaches \£i s \/ks, the deep impu- 
rities are all occupied but the rest of the population is shared 
by shallow localized and conduction band states. Initially, the 
conduction band electrons have a longer spin lifetime because 
impurity scattering is frequent at low temperatures so the DP 
mechanism that relaxes them is not very effective. However, 
the DP mechanism increases rapidly as T increases and the 
Z s curve turns around. At T > |e; s |/^b, the shallow impu- 
rity level is empty and the relaxation is dominated by the DP 

mechanism in the conduction band: 1 /i c = l/z^ p (T). 

At this point it is necessary to point out why only the linear- 
in-T DP mechanism is needed to explain the observed con- 
duction spin relaxation. The other two viable candidates (cu- 
bic DP and EY) for relaxation are much too weak to explain 
the observed relaxation times in ZnO. We use the calculated 
values for a^ P (th) and a^ P (th) in the previous section to ob- 
tain the relative relaxation efficiencies between the linear and 
cubic DP mechanism terms: 



l/x (1 



DP 



a { o' P (th) _1.73xl0 5 K 2 



V l DP 



a { ol(th)T 2 



r- 



(23) 



which demonstrates that the efficiency of the cubic-in-T term 
does not become comparable to the linear-in-T term at tem- 
peratures below 416 K which is far above the temperature 
range investigated here. In fact the cubic-in-T term does not 
even reach one-tenth the efficiency of the linear-in-r term in 
the temperature range investigated here. For this reason we 
can confidently ignore the cubic-in-T DP mechanism term in 
our fit. The crystal structure of ZnO therefore makes its spin 



relaxation qualitatively different from spin relaxation in bulk 
n-GaAs. We also compare the efficiencies of the DP and EY 
mechanisms: 



a^ P {th)x 2 p {T) 



1 / L DP _ 

1/tey a EY {th)T 



7.5 x 10 15 T 2 (r) Kns~ 



(24) 



Even if the momentum relaxation time taken to be unrealisti- 
cally low, say 1 fs, the DP mechanism is still nearly two or- 
ders of magnitude more efficient at relaxing spins than the EY 
mechanism in the temperature range studied here. Due to the 
drastic qualitative and quantitative differences between relax- 
ation mechanisms, we have unequivocally determined the rel- 
evant conduction electron spin relaxation mechanism in ZnO. 




0.01 



T(K) 



FIG. 1: Plot of l/z s vs. temperature. Points are experiment of Ref. 
4. Dashed curve: [«/. s /(n c +«; s )](1/Xdm)- Dotted curve: [n c /(n c + 
ni s )]{\/xop)- Solid curve: total l/x s . «,„„, t = 6.0 x 10 14 cm~ 3 



^imp.d 



5.0 x 10 17 cirr 3 ,£/ s 



L imp : s 

-23 meV, and 8/j - 



-360 meV. 



The fit of theory to the experimental data is clearly very 
good. We found that no reasonable fit was possible us- 
ing only a single impurity level, though this worked very 
well for GaAs 6 , so we used two levels. A good fit by this 
method was possible by adjusting the coefficients a^ P (exp) 
and ocdm (<?.*£>), and the binding energies £i s ,Eid and concen- 
trations n, m „ s ,n{ m p d of the two donors, subject to the con- 
straint that the room temperature carrier density should equal 
the measured 4 value of 1.26 x 10 15 cm~ 3 . Qualitatively, one 
finds that n imp j > n imPiS and |e M | > \e ls \ to get the right or- 
der of magnitude of the relaxation at low T. Physically, the 
deep impurity spins are important because they relax the shal- 
low impurity spins by the DM mechanism, and the strength 
of the low T relaxation implies that the deep impurities must 
be quite numerous. Quantitatively, a least squares fit to 

the data of Ref. (0) yields a^exp) = 134.5 Kr 1 ns~ 2 , 
a DM (exp)n imp j = 0.06 ns" 1 , |e w | = 360 meV, \e h \ = 23 



meV, and n 

(i) 



imp,s 



6.0 x 10 14 cm- 3 . 



a DP (exp) is about four times larger than the theoretical 

value of a^ P (th) given above, possibly due to strain effects. 
We also note that the values of z p that we used were taken 
from a different sample. 
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If we take ni mp ^ to be near the highest values measured 
for the deep donor (see below) then aoM( ex p) = 12 x 10~ 20 
cm 3 ns~' is about one order of magnitude larger than the the- 
oretical estimate UoM{th) given above. In view of the very 
poor understanding of the impurity wavefunctions, and the ex- 
ponential dependence of (Xdm on the overlaps, this is perhaps 
not too disturbing. 

The presence of a shallow donor and a very deep donor has 
been seen in hydrothermally grown ZnO samples of the type 
investigated here. 9 ' 31 Donor concentrations up to nearly 5.0 x 
10 17 cm~ 3 ( riimp d) have been measured for donors 330 — 360 
meV (|£m|) deep . 9,31 i 32 Donors as shallow as 13 — 51 meV 
(|e/j|) have been measured 8 at lower concentrations ~ 5.0 x 
10 14 cm~ 3 (n, mp ,. s ). Comparison with our values indicates that 
the parameters extracted from the fit are very reasonable for 
this material. 

From this analysis, we predict that in ZnO samples with 
fewer deep impurities, the relaxation time at low temperatures 
can be increased. As the impurities of ZnO vary greatly be- 
tween different growth techniques 33 , this prediction could be 
tested by further optical orientation experiments on different 
samples. 

VII. CONCLUSIONS 

We have found that % s in bulk ZnO can be understood by 
invoking previously known spin relaxation mechanisms. The 
dominant mechanisms in the material turn out to be the DP 
(scattering) relaxation of the conduction electron spins for 
T > 50 K and the DM (anisotropic exchange) mechanism for 
the localized spins for T < 50 K. In addition, it is very impor- 



tant to include the cross-relaxation between localized and con- 
duction states previously proposed for GaAs. These physical 
ingredients explain quantitatively the relatively fast relaxation 
at low temperatures as being due mainly to the DM mecha- 
nism which in turn depends on having both deep and shallow 
impurity states. At high temperatures, the conduction states 
are dominant, and the DP mechanism gives an excellent fit to 
the data. The combination explains the very surprising non- 
monotonic temperature dependence of x s . 

Finally, there are two aspects of the data in Ref. (3) that we 
have not addressed here: the applied magnetic field depen- 
dences on the spin relaxation and the spin relaxation observed 
in ZnO epilayers. We plan on addressing the former issue 
in a future publication. As for the latter issue, the epilayers 
are doped three to four orders of magnitude higher than in the 
bulk case. At such high dopings, spin glass effects become im- 
portant and localized donor states coalesce to produce donor 
bands; we do not expect our theory to be applicable in such a 
regime. 

The theory has now been sufficiently developed that optical 
orientation experiments can actually serve as a characteriza- 
tion tool for doped semiconductors, giving information about 
the binding energies and concentrations of the electrically ac- 
tive impurities in «-type materials. 
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